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Abstract 

A  closure  depth  is  established  from  bedform  stress  relations  that 
determines  the  seaward  extent  of  equilibrium  profiles  on  beaches.  The 
closure  depth  becomes  a  function  of  sand  size  and  incident  wave  height  that 
provides  the  matching  conditions  for  the  conjoined  shorerise  and  bar-berm 
portions  of  equilibrium  profiles  defined  by  the  elliptic  cycloid  solution  of 
Jenkins  and  Inman  (in  press).  These  relations  have  been  coded  in  the 
VORTEX  mine  burial  model  and  field  tested  (Jenkins  et  al.,  in  revision). 
Additional  findings  from  this  research  are  described  in  our  2002,  2003  and 
2004  Annual  Reports  [Grant  Number  N000 14-02- 1-0232, 
http://www.onr.navy.mil]. 

Introduction 

Closure  depth  for  seasonal  beach  profiles  is  usually  taken  as  the  depth  at 
which  the  envelope  of  profiles  converge  or  close  e.g.,  where  the  depth 
change  vs  depth  decreases  to  a  common  background  error  (Figure  1)  (e.g., 
Kraus  and  Harikai,  1983;  Inman  et  al.,  1993).  When  observations  are 
limited  to  comparison  of  two  or  three  surveys,  the  closure  depth  becomes 
that  at  which  the  survey  lines  converge  with  depth,  a  point  of  some 
uncertainty  (e.g.,  Shepard  and  Inman,  1951;  Nordstrom  and  Inman,  1973; 
Birkemeier,  1985).  Here  we  define  closure  depth  hc  as  the  maximum  depth 
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at  which  seasonal  changes  in  beach  profile  are  measurable  by  field  survey, 
most  commonly  using  fathometers  (Inman  and  Bagnold,  1963). 

Because  of  its  importance  to  the  understanding  of  seasonal  changes  in 
beach  profiles,  closure  depth  in  various  terminologies  has  a  much  studied 
history  associated  with  measurements  on  ocean  beaches.  Shepard  and  La 
Fond  (1940)  showed  that  seasonal  beach  changes  began  in  depths  greater 
than  5  m  and  extended  through  the  surf  zone  to  the  beach  face,  with  the 
maximum  changes  occurring  between  depths  of  1.5  m  and  3.0  m  in  the 
vicinity  of  the  breakpoint-bar.  Studies  of  profile  changes  show  similar 
winter/summer/winter  cyclic  changes  that  are  repeated  each  year  on 
perennial  beaches,  i.e.,  those  with  sufficient  sand  to  form  equilibrium 
summer  and  winter  beach  profiles  (e.g.,  Inman  et  al.,  1993;  Wright,  1995; 
Cowell  et  al.,  1999).  In  general  sand  from  the  deeper  winter  breakpoint-bar 
is  transported  onshore  during  small  summer  waves,  building  a  steeper  beach 
face  with  a  wide  summer  berm.  Higher  waves  of  winter  move  sand  offshore 
eroding  the  summer  beach  face  and  forming  a  gently  sloping  winter  beach 
out  to  a  deeper  breakpoint-bar  (Shepard  and  La  Fond,  1940;  Shepard,  1950). 

The  asymptotic  nature  of  the  closure  of  seasonal  beach  profiles  seaward 
of  the  surf  zone  makes  an  accurate  determination  of  closure  depth  difficult. 
Therefore,  early  on,  a  three  year  study  was  undertaken  adjacent  to  the 
Scripps  Institution  of  Oceanography  at  La  Jolla  California  with  the  dual 
objectives  of  obtaining  accurate  measurements  of  sand  level  change  in  the 
zone  of  profile  convergence,  and  of  evaluating  the  accuracy  of  fathometer 
surveys  of  sand  level  change.  A  series  of  five  bottom  stations  were 
established  in  water  depths  of  2.7  m,  6.3  m,  10.0  m,  16.7  m  and  22.2  m 
below  MSL,  each  consisting  of  an  array  of  six  reference  rods  forced  into  the 
bottom,  with  the  protruding  portions  measured  by  divers  (Inman  and 
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Rusnak,  1956).  The  two  shallow  stations  were  lost  during  seasonal  changes 
that  alternately  buried  the  2.7  m  station  while  the  6.3  m  station  rods  were 
lost  during  a  winter  cut  in  sand  level  that  exceeded  61  cm.  The  three  deeper 
reference  rod  stations  were  visited  and  measured  by  divers  and  the  stations 
surveyed  by  fathometer  on  a  fortnightly  to  monthly  basis  during  the  three 
year  period  of  1953-1955  (Table  1).  The  reference  rod  study  had  a  standard 
error  in  measurement  of  1.5  cm,  and  showed  that  seasonal  sand  level 
changes  occurred  out  to  depths  of  10  m  where  the  range  in  level  was  about 
10  cm.  The  two  deeper  stations  had  changes  of  5  cm,  which  were  not 
seasonal  in  nature,  but  due  to  shorter  period,  unspecified  phenomena. 

Inman  and  Bagnold  (1963)  in  their  analysis  of  seasonal  sand  level 
changes  in  the  energy  profiles  of  equilibrium  beaches,  also  show  the 
maximum  closure  depth  off  Scripps  Institution  of  Oceanography  to  be  about 
10  m.  They  noted  that  although  wave-induced  sand  motion  was  observed  to 
occur  out  to  52  m  depth  (Inman,  1957),  sand  motion  per  se  may  be 
associated  with  phenomena  other  than  the  closure  depth  of  seasonal  beach 
profiles. 

More  recently  Hallermeier  (1978,  1981)  derived  a  relation  for  closure 
depth,  by  assuming  that  closure  depth  is  related  to  sediment  suspension 
energetics  given  by  a  critical  value  of  the  Froude  number, 

hc=2.2Wa-6.S5(Hi/gT1)  (1) 

where  Hss  is  the  nearshore  storm  wave  height  that  is  exceeded  only  12  hours 
each  year  and  T  is  the  associated  wave  period.  Hallermeier  also  defined  a 
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deeper  depth  h\  as  the  limiting  depth  for  significant  on-offshore  transport  of 
sand  by  waves  throughout  a  typical  year, 


hl=(H,s-03ajT(g/5mD)u2  (2) 

where  ass  is  the  standard  deviation  of  Hss  and  D  is  the  characteristic  grain 
size,  usually  from  a  depth  of  about  1.5  hc. 

Birkemeier  (1985)  suggested  different  values  of  the  constants  in  equation 
(1)  and  found  that  the  simple  relation  hc=  1 .57  Hss  provided  a  reasonable  fit 
to  his  profile  measurements  at  Duck,  North  Carolina.  Cowell  et  al.  (1999) 
review  the  Hallermeier  relation  for  closure  depth  hc  and  limiting  transport 
depth  h\  and  extends  the  previous  data  worldwide  to  include  Australia.  Their 
calculations  indicate  that  hc  ranges  from  5  m  (Point  Mugu  California)  to 
12  m  (SE  Australia),  while  h\  ranges  from  13  m  (Netherlands)  to  53  m  (La 
Jolla  California).  They  conclude  that  discrepancies  in  data  and  calculation 
procedures  make  it  “pointless  to  quibble  over  accuracy  of  prediction”  in  hz 
and  hx.  In  the  context  of  planning  for  beach  nourishment,  Dean  (2002) 

observes  that  “although  closure  depth . is  more  of  a  concept  than  a  reality, 

it  does  provide  an  essential  basis  for  calculating  equilibrated... beach 
widths.” 

The  Hallermeier  relations  did  not  fit  our  modeling  needs  (Inman  et  al., 
2005a, b),  leading  us  to  seek  a  closure  depth  criterion  that  relates  bedform 
type,  sand  size  and  wave  height  as  outlined  below. 
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Bedforms  Associated  with  Equilibrium  Beach  Profdes 

Field  studies  of  bedform  and  type  of  sand  motion  indicate  that  their 
relation  can  be  expressed  in  terms  of  the  wave  form  of  the  Shield's  stress 
ratio,  9  -fl{ps  -  p)gD,  where  the  stress  t  is  taken  as  pum,  ps  and  p  are 

the  densities  of  solid  grains  and  water  respectively,  g  is  acceleration  of 
gravity,  D  is  median  grain  diameter,  and  um  is  the  near  bottom  amplitude  of 
the  orbital  velocity  (e.g.,  Dingier  and  Inman,  1976;  Conley  and  Inman, 

1992).  The  onset  (threshold)  of  grain  motion  6*t  =  (1.7  sec-1  )T  is  frequency 

dependent,  ranging  in  value  from  8.5  to  34  for  wave  periods  T  of  5  s  and 
20  s  respectively, 


Vortex  ripples . 0t  <  0  <  40 

Transition  ripples . 40  <6  <  240 

Flat  bed  (carpet  flow) . 9C  >  240 


The  above  limits  for  ripples  morphology  are  graphed  in  Figure  2,  where  the 
critical  data  from  significant  waves  Hs  are  shown  in  red  to  emphasize  the 
governing  importance  of  higher  waves  in  ripple  morphology. 

Care  is  to  be  exercised  in  the  use  of  Hr ms  vs  Hs  in  estimating  bedforms 
from  the  wave  form  of  the  Shields  parameter,  where  the  wave  stress  is 

assumed  to  have  the  form  f  =  pu2m,  without  the  usual  ambiguity  of  a  drag 
coefficient.  The  mechanics  of  grain  onset  <9t  and  ripple  formation  and 

disappearance  were  originally  based  on  the  dynamics  of  steady  state 
hydraulic  flows  in  laboratories.  Dingier  and  Inman  (1976)  show  that  it  is  the 
higher  groups  of  ocean  waves,  represented  by  Hs  that  are  analogous  to  the 
steady  state  laboratory  case.  In  the  thermodynamic  analysis  of  Jenkins  and 
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Inman  (in  press)  where  Hrms  is  applied  to  energetics,  it  is  the  occasional 
higher  waves  represented  by  Hs  that  form  and  modify  the  ripples.  Since 
Hs  =  a/2 Hrms  and  G  ~  u2m,  there  is  a  factor  of  2  difference  in  the  values  of 

G  calculated  from  the  two  wave  statistics,  but  only  that  from  Hs  gives  a 
meaningful  value  of  G .  In  summary,  for  any  given  ocean  wave  condition, 
ripple  morphology  is  a  function  of  the  significant  waves  Hs,  while  the  rates 
of  formation  and  sediment  transport  are  calculated  from  Hms. 


Vortex  Ripples 

The  vortex  ripple  regime  Gt  <G  <  40,  with  a  mid  range  of  about 

G  =  25,  implies  the  existence  of  a  rippled  sand  bottom,  where  the  groups  of 
high  waves  cause  stresses  above  the  threshold  of  motion  6{,  but  usually 
below  the  change  to  the  transition  ripple  regime.  Stresses  due  to  lower 
waves,  between  groups,  are  usually  below  Gt.  The  sand  would  be  motionless 
until  groups  of  larger  waves  suspend  sand  as  vortex  plumes  over  the  rippled 
bed.  In  vortex  ripples  there  is  no  instantaneous  granular  exchange  from 
ripple  to  ripple  unless  the  bottom  slope  is  sufficient  to  induce  a  gravitational 
component  of  sand  motion  over  a  horizontal  distance  of  about  A/2  or  greater. 
However  there  could  be  a  net  ripple  motion  associated  with  differences  in 
on-offshore  orbital  velocities  and/or  bed  slope  (e.g.,  Inman  and  Bowen, 
1962).  Generally,  this  would  be  a  small  slope  condition  of  zero  net 
transport,  seaward  of  the  closure  depth. 
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Transition  Ripples 

The  transition  ripple  regime  40  <6  <  240,  with  a  mid  range  of  0  =  140, 
indicates  a  sand  bottom  of  active  ripples  with  changing  ripple  steepness. 
Smaller  waves  between  groups  of  high  waves  form  vortex  ripples 
(T/77  ~  6)  with  distinct  vortex  plumes  above  the  ripple  crests,  while  groups 
of  higher  waves  reduce  the  ripple  height,  forming  transition  ripples 
(T/77  ~  12).  The  vortices  associated  with  transition  ripples  are  intense  and 
remain  closer  to  the  lee  of  the  ripple  crest.  Transition  ripples  occur  in  the 
deeper,  seaward  portion  of  the  shorerise  and  may  be  a  common  feature  in  the 
vicinity  of  closure  depth  for  equilibrium  beaches. 

Carpet  Flow  (flat-bed) 

Carpet  flow  is  the  most  active  and  intense  regime  of  sand  motion  under 
waves  and  occurs  for  Shield’s  numbers  above  the  critical  value  6C  =  240. 

The  bed  is  flat  and  wave  stresses  move  sand  as  sheets  or  carpets  of  moving 
sand  that  include  visually  distinct  forms  referred  to  as  streaking ,  roiling  and 
pluming  (Conley  and  Inman,  1992).  Carpet  flow  commonly  occurs 
throughout  the  inner  portions  of  the  shorerise  and  in  the  vicinity  of  the 
breakpoint  of  waves. 

Closure  Depth 

While  it  may  be  reasonable  to  apply  (1)  and  (2)  or  its  simpler  form  after 
Birkemeier  (1985),  comparisons  with  the  Inman  et  al.  (1993)  beach  profile 
data  set  show  that  these  relations  tend  to  underestimate  closure  depth.  We 
propose  an  alternative  closure  depth  relation.  This  relation  is  based  on  two 
premises:  closure  depth  is  the  seaward  limit  of  non-zero  net  transport  in  the 
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cross-shore  direction;  and,  closure  depth  is  a  vortex  ripple  regime  in  which 
no  net  granular  exchange  occurs  from  ripple  to  ripple.  Inman  (1957)  gives 
observations  of  stationary  vortex  ripples  in  the  field  and  Dingier  and  Inman 
(1976)  establish  a  parametric  relationship  between  dimensions  of  stationary 
vortex  ripples  and  the  Shield's  parameter  0  in  the  range  3  <  0  <40.  Using 
the  inverse  of  that  parametric  relation  to  solve  for  the  depth  gives, 


KeHa 


sinh  kh 


'd* 

c  vAy 


(3) 


where  Ke  and  yj  are  nondimensional  empirical  parameters,  D2  is  the  shorerise 
median  grain  size;  and  D0  is  a  reference  grain  size.  With  Ke  ~  2.0,  yj  ~  0.33 
and  D0  — 1 00  fj, m,  the  empirical  closure  depths  reported  in  Inman  et  al. 

(1993)  are  reproduced  by  (3).  Figure  3  gives  a  contour  plot  calculated  from 
(3)  showing  the  rates  at  which  closure  depth  increases  with  increasing  wave 
height  and  decreasing  grain  size.  Because  of  the  wave  number  dependence, 
closure  depth  also  increases  with  increasing  wave  period.  Figure  3  is  based 
on  T  =  15  sec.,  typical  of  storm  induced  waves  on  exposed  high-energy 
coastlines. 


Cycloid  Theory  and  Evaluation  of  Boundary  Conditions 

Understanding  of  the  foregoing  measurements  of  closure  depth  can  now 
be  placed  in  the  more  rigorous  framework  of  theory.  Jenkins  and  Inman  (in 
press)  show  that  shorerise  and  bar-berm  equilibria  were  found  to  have  an 
exact  general  solution  belonging  to  the  class  of  elliptic  cycloids  that  admit  to 
a  simple  analytic  form  h  =  Axm,  where  h  is  depth  taken  positive  downward 
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from  mean  sea  level  and  *  is  the  cross-shore  component,  positive  in  the 
offshore  direction.  The  elliptic  cycloid  formulation  gives  an  exact  general 
solution  that  tells  us  the  nature  of  the  equilibrium  beach  profile  curve  and 
provides  valuable  guidance  to  understanding  and  matching  of  conjoined 
curves  to  empirical  boundary  conditions.  The  eccentricity  e  of  the  ellipse  is 
the  most  fundamental  variable  in  the  solution.  It  is  governed  by  the  stress 
nonlinearity  of  the  waves  in  the  shorezone,  which  in  turn  is  a  function  of  the 
height  and  period  of  waves,  the  morphology  of  the  inner  shelf,  and  the  size 
of  sediment. 

The  equi  librium  profile  consists  of  two  conjoined  sections,  joined  at  the 
breakpoint  of  the  waves  (Figure  1).  The  bar-berm  inner  curve  extends  from 
the  berm  crest  above  the  beach  face,  through  the  surf  zone  to  the  breakpoint- 
bar;  and  the  shorerise  outer  curve  extends  from  the  wave  breakpoint  seaward 
to  the  closure  depth.  The  elliptic  cycloid  theory  shows  that  both  the  bar- 
berm  and  the  shorerise  have  similar  curves,  each  defined  by  the  trace  of  the 
rotating  ellipse  (Figure  4).  The  two  curves  have  the  same  cycloidal  form  but 
may  differ  in  detail  depending  upon  the  eccentricity  e  of  the  ellipse  and  the 
boundary  conditions. 

Before  proceeding  with  the  analysis  of  closure  depth  it  is  of  interest  to 
consider  the  larger  family  of  which  the  elliptic  cycloid  is  a  special  case.  By 
analogy  to  the  cycloid  traced  by  a  point  on  the  surface  of  a  rolling  circle, 
which  is  a  special  case  of  a  trochoid,  elliptic  cycloids  are  the  special  case 
where  the  cycloid  is  traced  by  a  point  on  a  rolling  ellipse.  This  family  of 
curves  also  includes  elliptic  trochoids  traced  by  points  that  extend  outside  of 
the  ellipse  on  a  semimajor  or  minor  axis  (prolate  elliptic  cycloids),  as  well  as 
on  an  axis  within  the  ellipse  (curtate  elliptic  cycloids).  Out  present,  first, 
treatment  of  theoretical  equilibrium  beach  profiles  addresses  only  the 
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common  elliptic  cycloid,  where  the  point  tracing  the  cycloidal  curve  is  fixed 
on  the  circumference  at  the  end  of  the  semiminor  axis  (Figure  4).  It  is  to  be 
noted  that  the  slope  at  the  seaward  (flat)  end  of  the  profile  (6>=  tt)  is  always 
zero  for  all  members  of  the  elliptic  trochoid  family  of  curves.  However  the 
profile  slope  at  the  landward  end  (9=0)  varies  with  type  of  elliptic  cycloid 
(prolate/common/curtate).  In  the  future,  this  refinement  could  be  addressed 
in  conjunction  with  the  slopes  of  the  upper  beach  face  in  the  bar-berm  and 
the  seaward  slope  from  the  breakpoint  of  the  shorerise. 

In  the  bar-berm  curve,  the  zero  slope  in  the  elliptic  cycloid  occurs  at  the 
breakpoint-bar,  at  a  depth  of  hb.  However,  the  zero  slope  of  the  shorerise 
curve  at  closure  depth  suggests  that  this  is  the  depth  of  water  necessary  for 
an  equilibrium  profile  to  exist,  here  called  the  theoretical  closure  depth  hc. 
Along  ocean  beaches  with  seaward  sloping  shelves,  the  necessary  depth 
always  exists  but  usually  with  a  measurable  slope  f>  0. 
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Table  1.  Summary  of  sand  level  changes  on  the  beach  and  shelf  at  La  Jolla,  California. 
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Figure  1.  Envelope  of  variability  of  12  measured  beach  profiles 
(1950-1987)  at  Oceanside,  CA  (gray)  compared  with  the  ensemble 
of  type-a  elliptic  cycloid  solutions  (colored)  for  selected  incident 
wave  heights  and  periods  with  shorerise  sand  size  D2  =  100/zm  and 
bar-berm  sand  size  D{  =  200  /urn.  Colored  lines  break  at  conjoined 
point  and  end  at  closure  depth  (from  Jenkins  and  Inman,  in  press). 
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Figure  2.  Ripple  steepnes  r| / A,  vs  relative  stress  0  for  steady  laboratory  (O,  A)  and  ocean  waves  (□  ,  ,  :  ), 

the  latter  computed  from  significant  wave  height  Hs.  Data  trends  shown  as  dashed  lines;  short  vertical  lines 

give  onset  of  grain  motion  0t  for  selected  wave  periods  T  (after  Dingier  and  Inman,  1976).  04.k9e-4a 
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Figure  3.  Closure  depth  hc  from  (3)  showing  dependence  on  shorerise  grain 

size  D2  and  incident  wave  height  for  waves  of  1 5  sec  period 

(i//  =  0.33,  Ke  -  2.0 ,D0  =  100//m)  (from  Jenkins  and  Inman,  in  press). 
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Figure  4.  Elliptic  cycloid  solution  (red  line)  for  the  equilibrium  beach  profile  traced 
by  a  point  on  a  rolling  ellipse  with  semimajor  and  semiminor  axis  a,  b,  eccentricity 

e-^fl- ( b 2  / a 2 ),  polar  coordinate  vector  r,  angle  of  rotation  6,  and  horizontal  length 
scale  stretching  factor  equal  1.0  (from  Jenkins  and  Inman,  in  press). 
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